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Abstract. 

In the limit d ^ oo the role of pressure gradients and that of the incompressibility 
constraint decreases, thus blurring the difference between transverse and longitudinal velocity 
correlation functions. Using Polyakov's expression for the dissipation anomaly the closed 
equation for the probability density function is obtained. This model for the dissipation 
terms is the only one satisfying both equations of motion and a set of dynamical constraints. 
The resulting equations show that when d ^ oo, the predictions of Kolmogorov theory are 
exact. It is also shown that the 0{l/d) pressure effects, producing the regularization of the 
equations of motion for the PDF are responsible for the distinction between the Navier-Stokes 
and Burgers dynamics. 
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The problem of d- dimensional Navier-Stokes turbulence was first considered by Frisch, 
Fournier and Rose [1] who hoped to develop a renormalized perturbation expansion with 
1/d as a small parameter. This attempt failed because the 1/d-factors appearing in the 
Wyld diagrammatic expansion due to the angular integrations are cancelled by the 0{d)- 
multipliers, resulting from the summation over the d components of the velocity field. In 
this paper, using the nonperturbative approach, we will revisit the problem of d-dimensional 
turbulence 

The equations of motion are (density p = 1): 

dtVi + VjdjVi = -diP + vV'^Vi + fi (1) 

and 

diVi = (2) 
where f is a forcing function mimicking the large-scale turbulence production mechanism 



and in a statistically steady state the mean pumping rate P — f • \ — £ — 0{1) where S is 
the mean dissipation rate. In what follows we will also deal with the local dissipation rate 



so that E — S{x). The mean dissipation rate is independent on the space dimensionality d. 

Among other results, Kolmogorov theory produced a clear distinction between longitudi- 
nal and transverse structure functions: Consider two points x and x' and define r = x — x'. 
Assuming that the x-axis is paralel to the displacement vector r, one can find that in the 
inertial range [2], [3]: 

dy+^S^ = -^£ (3) 



giving 



12 



^3 = {/^uf ^ {u{x') - u{x)f ^ -^^j^f (4) 



and 



SI = {Avf = {v{x') - v{x))^ « (5) 
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where u and v are the components of velocity field paralel and perpendicular to the a:-axis 
(vector r). These two relations show that the energy flux in the d-dimensional turbulence 
is dominated by the longitudinal velocity fluctuations only and that the remaining d — 1 
transverse components do not directly participate in the energy transfer. It is easy to show 
that the probability density P{Au,r) = P(— Aw, — r), while P{Av,r) = P{—Av,r). Thus, 
the time-reversal symmetries of transverse and longitudinal velocity fluctuations are very 
different. This difference is lost in the Fourier-transforms of the relations (4) and (5), both 
giving S3{k) — Sl{k) = for A; 7^ 0. That is why the Wyld expansion, deahng with all 
components of velocity fleld on an equal footing, is unable to make a crucial distinction 
between longitudinal and transverse velocity fluctuations. The relations (4), (5) also show 
that in the limit d — > 00 the asymmetry of the probability density function, responsible for 
non-zero values of the odd-order moments decreases. 
We consider the A'"-point generating function: 

Z =< e^*-"("') > (6) 

where the vectors Xi define the positions of the points denoted 1 < i < N. Using the 
incompressibility condition, the equation for Z can be written: 

with 

If = Y.<^y f(xj)e^''^(^^) > (8) 
3 

7p = -^A,<e^'"(^')^^> (9) 
■ dxj 

and 

D = uJ2>^3< e^»"(^>)^^?^if^ > (10) 
j "^3 

In what follows we will be mainly interested in the probability density function of the 
two-point velocity differences which is ontained from (7)-(10), setting Ai -|- A2 = (sse Refs. 
[4]- [6]), so that 

Z =< exp{X ■ U) > (11) 



where 



U = u(x') - u(x) = Au (12) 

The moments of the two-point velocity differences which in homogeneous and isotropic turbu- 
lence can depend only on the absolute values of two vectors (velocity difference v(x') — v(x) 
and displacement r = x' — x) and the angle 9 between them with 9 — 7r/2 and ^ = corre- 
sponding to transverse and longitudinal structure functions, respectively. It is easy to show 
[2]- [3] that the general form of the second-order structure function in the inertial range is: 

S2{r,9) = ^^D,L{r){l - ^^cos'{9)) (13) 

with DLL{r) —< {u{x) — u{x + r)y >. More involved relation can be written for the fourth- 
order moment: 

S^{r, 9) = DLLLL{r)cos\9) - 3DLLNN{r)sin'{29) + Di,^N^{r)sin\9) (14) 

where D^lnn —< {v{x) — v{x + r)y{u{x) — u{x + r))^ > and v and u are the components 
of the velocity field perpendicular and parallel to the x-axis, respectively. In general, in the 
Uimit cos{9) = s ^ ±1, corresponding to the moments of the longitudinal velocity differences 
Sn{r, s) Sn{r)cos"'{9). This means that in this limit Z(A, r, s) — >• Z{Xs, r) = Z{Xx, r). The 
generating function can depend only on three variables: 



A ■ r I — 
Vi ^r; 7/2 = = Xcos{9); 773 = yA^ - 77I; 



The equations (7)- (10) become after some manipulations: 

Zt + [dA + ^9,2 + -d,As + -dmdm + ^^^^^,3 - ^^^jZ = I, + D (15) 
where 

Ip = K< {d2,iP{2) - ai,,p(l))e^ u > (16) 

and 

D = Ky< {dl,v,{2) - dl,v,{l)e^-'' > (17) 
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where, to simplify notation we set di^a = and v{i) = v(xi). Interested in the hmit 
773 —> 0, we will seek a solution in a form: 



'Hi 

The equation for the generating function becomes (the subscript o is omitted hereafter): 



Z ^ Z,(r)„m)K{{^n (18) 

V2 



[dm9m + 9r,, + —d^, + ]Z ^Ip + D 19 

r ' r]2 r]2 r 

where a = K'{0). When rj^ —>■ 0, one can introduce a more general relation Z a; Zo{r)i, 772) + 
771^1(771,772) and derive an equation for Zi [7]. This, however, does not give any advantage 
over the expression (18). The turbulence production contributions, which in the case of the 
large-scale forcing function rapidly varying in time can be written as: 

If = ^^2^1(1 + 2s') ^ ^iPvlZ (20) 

will be discussed later. The form of Z, introduced above, is only an approximation, valid in 
the d ^ 00 limit. One can rewrite the dissipation term 



D = -M < (£(1) + £(2))e''^^ >+s= - ^Z - ^rjl < [6{1) + 5(2)]e''^^ > +s (21) 

where the dissipation fluctuations are defined as 5{x) = 8{x)—£. The "surface contribution" 
s: 



in the large Reynolds number limit 1/ — > 0. Substituting this into (19) gives after Fourier 
transformation: 



-duUdrP-- — %uUP+2adrP+2a- — -P = ^dfjP{U)+-^dljK{U)P{U)+d^UU)P{U) 
r r d d 

(22) 
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where 

k{U)^<{Si + S2)\U > (23) 

and 

ip=<d2,M'^)-9i,M^)\U> (24) 

are conditional expectation values of Si + 82 and 92,a;p(2) — di^^Pi^) for a fixed value of the 
longitudinal velocity difference U = Ux, respectively. It is clear that 

/cxD rco 
K{U)P{U)dU = / ip{U)P{U)dU =< Si + S2 >=< ip >= (25) 
-00 J —00 

Multiplying (22) by gives: 

^ a,r('='-i)-53 = -l^ + i (26) 



where 

/oo 
ip{U)UP{U)dU (27) 
-00 

a=i^ (28) 
1 + 2q;/3 ^ ' 

This result differs from exact Kolmogorov relation (3), involving the factors r''+^ instead 
of r'^~^ derived above. The error stems from the shape of the generating function (18), 
introduced above. In reality the transverse and longitudinal correlation functions are related 
by the incompressibility constraint 

drS2 + S2^Sl (29) 



d~ 1 

providing a boundary condition on a mixed derivatives of the Z-function at A = 0. This 
relation is not accounted for in the anzatz (18). One can see that when d — > 00, the 
role of the incompressibility constraint and, as a consequence, of the pressure gradients, 
decreases leading to the longitudinal and transverse correlations functions very close to each 
other. This is easily understood since a single constraint on the d — > 00 velocity components 
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cannot quantitatively change too much. It will become clear below that although, in general, 
the effect of the pressure contributions is 0(l/ci), they cannot be neglected since they are 
responsible for the regularization of the solution, thus producing a major difference between 
the Navier-Stokes and Burgers dynamics. Comparing (26) with (3) we conclude that when 
d is large a = 1, i = and that the the solution in the form (18) is a good approximation. 
In general we know nothing about the function k{U) which is a part of the dissipation term 
D. Luckily, it is irrelevant in the above calculation of 5*3. To obtain a closed equation for 
the generating function, the expression for D is needed. It has been shown by Polyakov [4] 



where we choose A — + c and 7 and c are 0{1). The relation (30) looks the same as 
the one appearing in the Polyakov theory of Burgers turbulence. This is not accidental 
since, barring the pressure conntributions, the equation of motion (1) is formally similar 
to the multidimensional Burgers equation. The expression (30 ) is a result of an operator 
product expansion, based on a point-splitting calculation, developed for the one-dimensional 
problem of Burgers turbulence [3]. It is easy to show that (30) (see Refs. [4], [5]) exactly 
satisfies the equation of motion (1) (without pressure) and is the only possible model, not 
violating scaling (see below) and both Galilean invariance and an obvious constraint U — 0. 
Substituting this into (22 ) with a = 0, we have after Fourier transformation: 



(30) 



duUdrP + 



{d-1) 



duUP = -AdrP - dl 



ip{U)P{U) 



(31) 



r 



Keeping only the 0(d)- contributions gives: 



- -duUP = ^ddrP 



(32) 



r 



Multiplying this by t/" and integrating gives: 



n-Sn = 'yddrS, 
r 



n 



(33) 



Seeking the solution in a scaling form: Sn oc r^" gives 
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Cn ^ n/3 (34) 

where 7 = 3 is chosen to satisfy the relation (1). Thus, the predictions of Kolmogorov theory 
are vahd in the hmit d ^ 00. One can see that in this case 

P(U,r)^\F{^) (35) 

f-i r 3 

which is consistent with the scahng assumption involved in derivation of the expression for 
the dissipation anomaly (30). This result needs some explanation. The d — > 00 limit must 
be understood in a following way: first we fix the moment number n and then drive d to 
infinity. It is clear from (31) that neglecting the pressure term the exponents 

= S (36) 
are singular at n = ^4 and violate all possible dynamic constraints. This happens due to the 
sign of the 0(1) contribution to the equation of motion (22). Thus, the role of the pressure 
terms is in modifying the sign of the "Burgers- like" first contribution to the left side of 
(22), thus producing the regularization, distinguishing Burgers dynamics from that of the 
Navier-Stokes. The possibilities of modelling the pressure contributions will be discussed 
in a future communication. It is interesting, that gradual disappearence of the anomalus 
scaling was observed in a numerical solution of a set of N coupled one-dimensional "shell 
models" : the deviations from Kolmogorov scaling decreased to close to zero with increase of 
TV [8]. 

Now, we can show that the pressure terms Ip — 0{1) and are small compared to the 0{d) 
contributions. The generating function in the limit d — > 00 can be written as: 

Z ^ U-^TaJ-^ (37) 

Assuming that at the integral scale r — L — 1 the PDF becomes close to the gaussian gives 
A2„ = 2-(2n-l)!!: 
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Z =< e^^ >;:^ e^^''-^ 

and 

V< e2^^ > < Z 

since X — ik with real /c. Using this relation the Schwartz inequality is: 

Ip < ^{X-V^^v^'^^ Z (38) 
The four- order correlation function VaVpVjVs ~< VaVj >< VfjVs > +.... where 

with = 0(1). It means that VaVaVp — — 0. The only sorce of the 0{d) 

contributions is V^rj = d which cannot appear in (38) due to incompressibility. That is why 
Ip < QZ, where Q is the operator, involving rji, r]2 and derivatives, all 0(1). 
To conclude: Kolmogorov scaling is an exact solution of the Navier- Stokes equations in 
the limit of large space dimensionality d ^ oo, provided the assumptions leading to a self- 
consistent conjecture for the dissipation anomaly, similar to the one derived by Polyakov [4] 
in the theory of Burgers turbulence, are valid. We do not know if this solution is the only 
possible one. The 0{l/d) pressure contributions to the Navier-Stokes equations, though 
unable to strongly influence the values of the scaling exponents, are crucial providing the 
regularization of the otherwise singular expression for the scaling exponents. These terms 
are responsible for a clear distinction between the Navier-Stokes and Burgers dynamics. 
The equation (7) can be written for the multidimensional Burgers equation if density fluc- 
tuations are taken into account [4], [5]. The fact that in the limit d oo, the pressure 
terms are unimportant make the Navier-Stokes and Burgers equations formaly seem the 
same. This not so. The steady state solution in the random- force-driven Burgers equation 
is impossible without accounting for the forcing term (20) which imposes a typical Burgers 
scahng behaviour U — 0{r) [4], [5]. This scaling is not always there: It has been shown in [9] 
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and [10] that, in principle, the scahng of velocity correlation functions in Burgers turbulence 
depends on the properties of the forcing function. In the case of the Navier-Stokes dynamics 
one can drop the forcing term and obtain the inertial range Kolmogorov scaling U — 0{r^). 
However, as was shown above, neglecting the pressure contributions leads to the singularity 
of some high-order moments and to the unphysical behaviour. Still, when d is large, due to 
formal similarity of the equations of motion, one cannot rule out Kolmogorov scaling as a 
short-time asymptotics of decaying Burgers turbulence. 
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